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Abstract 

The nature of the electromagnetic (EM) energy for general charge 
and current distributions is analyzed. We come to the conclusion that 
EM energy can reside only within charge-current distributions, and 
that EM energy density cannot exist in otherwise empty space. The 
form ^{p(f>+yA) is a suitable EM energy density, but (E-D + B'H)/87r 
is not. 

PACS numbers: 03.50.De, 41.20.-q 



1 Introduction 

James Clerk Maxwell, in his 1864 paper "A Dynamical Theory of the Elec- 
tromagnetic Field" [1], introduced the concept of electromagnetic (EM) field 
energy proportional to (E-D + B-H) that existed throughout space. For 
Maxwell this was a natural property of the aether which he conceived of as 
a material substance permeating space. He thought of the EM energy in the 
aether as being similar to elastic energy in any other material object. 

Although most physicists have long discarded the concept of an aether, 
at least as a material substance, the tendency has been to follow Maxwell's 
lead in the use of EM field energy. While most texts point out the alter- 
native of EM energy residing in the charge and current distributions in the 
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form of an energy density proportional to (p0 + j-A), they usually settle on 
(E-D + B-H) as the more useful form. Indeed, we have become comfortable 
with EM field energy permeating otherwise empty space. 

In this paper, we first analyze EM energy for specific examples of elec- 
trostatic (Sec. 2) and magnetostatic (Sees. 3 and 4) energy, and then gen- 
eral time dependent electromagnetic energy (Sec. 5). In each case we find 
that the EM energy must reside in matter, and reach the conclusion that 
(E-D + B-H)/87r cannot be considered an energy density. EM energy re- 
sides in charge and current densities. Only |(p0 + j-A) is a suitable EM 
energy density, with no electromagnetic energy in otherwise free space. 

2 Electrostatic Energy 

The electrostatic energy of a charge distribution p in a static electric potential 
00 is given by 




with the understanding that the external potential 0o is due to charges other 
than the distribution p. The subscript p0o is used here to emphasize that the 
energy given by this equation is due to a charge distribution in an external 
potential. This equation can be considered the definition of the electric 
potential as potential energy per unit charge. Note that the integral is taken 
over all space, but only regions with charge contribute to the integral. 

The electrostatic energy can also be given in terms of the electric fields 
by using Maxwell's equatioiil] for V-E' in Eq. ([T]): 




= i- /0o(V-E,)dV 

= ^/ [V-(0oE,)-E,-V0o]c^'r 

= i_ ^ dS-(0oE,) + ^ 1 E,-EorfV. (2) 

The surface integral is usually discarded by being taken over an infinite sphere 
^We do not include polarizable matter in this and the foUowing section. 
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for which it is presumed to vanish. This leads to the form 



t^EpEo = ^ / Ep-EorfV, (3) 

for the electrostatic energy purely in terms of the electric field Ep of the 
charge distribution and the external electric field Eq. (Note: We are using 
Gaussian units.) This integral is also over all space, however the regions of 
space that contribute to the integral are generally different for the integrals 
in Eqs. ([I]) and ([3]). 

The derivation of Eq. from Eq. ([T]) is so simple that they are generally 
considered to be equivalent formulas for the same electrostatic energy. How- 
ever we shall see that the two are not equivalent, and that Eq. ([3]) must be 
supplemented to give the correct energy for some cases. We first demonstrate 
this by the simple example of an electric dipole in a uniform electric field Eq. 
The potential for this electric field is 

00 = -r-Eo. (4) 

The energy Uptf,^ is given by 



U. 



P<f>0 



j p(j)od^r = - J pr-EorfV = -p-Eo. (5) 



We see that in a uniform electric field, the energy of a charge distribution 
depends only on its electric dipole moment p. For a point dipole at a position 
r in a spatially varying field, the energy is given by 

Up^ = -p-E(r), (6) 

For this case, the force on the point dipole will be the negative gradient of 
the energy: 

F = V(p-E). (7) 

Equation ([3]) using the electric fields gives a different answer for the en- 
ergy. For simplicity, we take the charge distribution to be a point dipole. 
Then the dipole electric field is given by 

E, = MEi^ _ 4.(p.r)M(r). (8) 
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When this is put into Eq. Q, the first term cancels in the angular integration, 
and the second term leads to 



'0 



J (p-f)(f-Eo)5(r)rfV 




(9) 



which is different than the energy Up^^. Of the two results, it is clear that 
— p-Eo is the correct energy, both from the definition of the potential and 
the fact that it leads to the correct form for the torque on the dipole. 

What went wrong with ?7epEo? We can see this by looking at the surface 
integral which is usually discarded, but does not always vanish. For some 
cases, such as the present one, the surface integral In Eq. ([2]) must be included 
to get the correct energy. For the electric dipole in a uniform field, the surface 
integral, taken over a sphere of radius r, is 



Including this surface contribution gives the correct electrostatic energy, 
in agreement with f/p^^ We see that the surface integral in Eq. ([2]) must be 
included except in those cases where it does vanish. Since the surface integral 
involves the potential 0, this means that the electrostatic energy cannot 
always be given solely in terms of the electric field. Note that the surface 
integral is independent of the radius of the sphere, which can be infinite. 
This means that even integrating E^ over "all space" does not exclude the 
need for the surface integral. 

This does not end the story, because there is an alternate derivation of 
Eq. ([3]) that seems unrelated to the potential. This derivation starts with 
the rate of energy put into matter by the electromagnetic field in a volume 
V, given by 



where j now represents all the currents and E is the electric field due to all 
sources. The standard textbook derivation, [2-4] using Maxwell's equations 




(10) 




(11) 
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and some vector identities, eventually leads to the result 



f 5^(e2 + B^)£r = -^i dS-(ExB). (12) 

dt Stt Jv 4:71 Js 

The surface integral on the right hand side is usually interpreted (because of 
its minus sign) as the rate in which electromagnetic energy enters the volume 
V through the surface S with the identification of the Poynting vector 

P=|-(ExB) (13) 

as the electromagnetic intensity (transmitted power per unit area). The first 
time derivative on the left hand side of Eq. (fT2l) is the rate of change of the 
energy of the matter inside the volume. Assuming conservation of energy, 
the second term is then deduced to be the rate of change of the energy of the 
electromagnetic fields inside the volume. For the purely electric case, this 
has been taken to imply that the electrostatic energy can be given by 

Ue = ^ f E'd'r. (14) 

The E in this integral is the total electric field, both external and that of 
the charge distribution p, so E^ = E^ + Eq + 2Ep-Eo. When the self energy 
contributions E^ and Eq are excluded from the integration, this gives the 
same result as Eq. ([3]). 

This will give the same wrong answer for the case of the electric dipole in 
a uniform electric field. What is missing now? To see this, let us first look 
at Eq. (ITTi) . the starting point for this second derivation,. We can consider 
the negative of the right hand side to be the rate of energy input to the 
electromagnetic field. If we start with a constant external field Eq and no 
initial charge distribution, the time integral of this rate will be the energy 
put into the field: 

U{t) = - [ dt I j-EorfV. (15) 

JO Jv 

Now, consider a rod of length L, originally uncharged. A current j in the rod 
in the direction of L will produce a dipole moment p(i(:) given by 

p(t) = g(t)L= f dt I ]d'r. (16) 
JO Jv 
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Using this expression for the electric dipole moment of the rod, we get for 
its energy in an external field Eg, 



U = - f'dt f j-Eoc/V = -p-Eo, 

Jo Jv 



(17) 



which is the correct answer. 

We see that the starting point for the derivation leading to Eq. (HM 
is correct. Where did the rest of the derivation go wrong? During the 
charging of the ends of the rod, there will be a current in the rod, which 
produces a magnetic field B. Thus there will be an ExB vector, resulting 
in electromagnetic energy passing through the outer surface of the volume. 
(This is if we want to use the E^ integral for the energy.) We now show 
that this energy equals two thirds of the energy put into the field by the j-E 
integral. 

We consider the rod of length L to be along the 2;-axis and use spherical 
coordinates. The magnetic field due to a current / in the rod is given by the 
law of Biot-Savart: 



B 



/0 



cr sin 9 



cos 9 



(r cos 6* — L) 



Vr^ + -2rL cos 9 



(18) 



This is the magnetic field produced by the current in the rod as its ends are 
being charged. According to Eq. (fT2l) . the electromagnetic power passing 
through a sphere of radius R during the charging of the rod in the presence 
of the electric field Eq is given by 



P 



^/dS.(E„xB)^^ 

-cR^ r ,^ IE, sin 9 



R^dQ 



f- (EoX0 



dn- 



Atx J cR sin 9 

-RIE, 



COS 9 



B 



{Rcos9 - L) 



du 



u 



VWTW^2RLcos9 
{Ru - L) 



y/R^ + L2 - 2RLu 



(19) 



We have left out a term in cos0 that integrates to zero. The remaining 
integral can be done by expanding in Legendre polynomials, giving 



-/Eq'L. 



(20) 
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This result is independent of the radius of the sphere. The energy that 
escapes beyond the sphere in producing the dipole moment (if we want to 
use the form for energy) is then 

[/escape = J Pdt = J lE^-Ldt = ~p-Eo (21) 

Thus two thirds of the energy escapes and only one third of the energy 
is represented by the volume integral of E^. This is consistent with what 
we found by doing that integral. The energy passing through the surface is 
accounted for by including the surface integral term in Eq. ([21). That is why 
the surface integral is necessary to give the correct electrostatic energy if the 
E^ form is used. 

We now consider the question of what combination of fields, potentials, 
or charges might be considered an electrostatic energy density. The electro- 
static energy in terms of the total potential and charge distribution is usually 
calculated by starting from having no charge distribution, and then building 
up the final charge distribution in infinitesimal steps. Using Eq. ([T]) for the 
incremental energy at each step, leads to the result 

Up4> = \j P<Pd'r: (22) 

The factor | arises from the relation 

5(p0) = p6(f) + (f)6p = 2(j)Sp, (23) 

which follows from p6(j) = (pSp in vacuum or a linear medium. The potential 
is now the potential due to all sources. 

Equation fl22|) gives the total electrostatic energy due to a charge distri- 
bution p if the integral is taken over all space. However, contributions to 
the integral come only from those regions where there is charge. If we inter- 
pret the integral over a finite region as giving the electrostatic energy in that 
region, we can define an energy density as 

"p./- = ^P(P- (24) 

With this definition, the electrostatic energy only exists in regions with a 
non-vanishing charge distribution. 
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Starting from Eq. (1221) . we can repeat the steps that led to Eq. ([2]) with 
the result 

= ^ = ^/'*«-<*^> + s/^'''''- 

for the electrostatic energy in terms of the total electric field. As before, 
the correct energy requires the inclusion of a surface integral involving the 
potential as well as E. 

What if the integral of is taken over a finite region? Can E^/Svr 
be considered an energy density? It is clear from Eq. (125|) that this is not 
the case unless the surface integral vanishes. Correctly including the surface 
integral results in the electrostatic energy being the same as that given by the 
integral over |p0. In any event, the first equality in Eq. fl25|) shows that there 
is no electrostatic energy in any region without charge. The fairly common 
belief that there can be electrostatic field energy in otherwise empty space is 
wrong. Our conclusion is that E^/Svr cannot be considered an electrostatic 
energy density. 



3 Magnet ost at ic Energy of a Current Distri- 
bution 

The energy of a current distribution in a external magnetic field Bq can be 
found by considering the hypothetical situation of turning on the magnetic 
field in the presence of the current, which is kept constant as the field is 
turned on. The constant current source that maintains the current provides 
the energy for this process. We start from Eq. (ITSl) for the energy input to 
the EM field: 

UiAoit) = - fdt' f j.E(t')rfV 

JO Jv 



dt' I Ud[Ait'))d\ 
Jv 

- I j-Ao(t)rfV. (26) 

c Jv 



'V 

In deriving this, we used the fact that the current j does not change with 
time, and have assumed the absence of an electrostatic potential. Equation 
fl26|) gives the energy in a magnetostatic field in terms of the current density 
and the vector potential. The integral is over all space, but only regions with 
current contribute to the integral. 
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We can introduce the magnetic field into Eq. fl26l) for tlie energy of a 
current distribution j in a static magnetic field Bq by using Maxwell's VxB 
equation: 

f/jAo = - / j-AorfV 

C JV 

= ^/,|V-(BjXA„) + Bj.(VxA„)l#r 

= s//S-(BiXA.) + ^|B,.BorfV, (27) 

where Bj is the magnetic field due to the current j. As was the case for 
the electrostatic energy, there is a surface term that is usually discarded in 
describing the magnetostatic energy in terms of B. Doing so leads to the 
expression 

^^^^'' = iX^^-^°^'" ^^^^ 
for the magnetostatic energy purely in terms of the magnetic fields. 

Equations ( l26l) and ( 127|) give the magnetostatic energy of a current dis- 
tribution in an external magnetic field. To derive the corresponding forms in 
terms of all currents and fields, we include an increase in the current in Eq. 
( !26|) utilizing the relation 

dt{yA) = -yidtA) + A-{dtS) = 2j-(9tA), (29) 

which holds in vacuum and within linear magnetic materials. This leads to 
the results 

U-At) = ^/^j-At/V (30) 

= — / dS-(BxA) + — / B^dV, (31) 

where the currents and fields are due to all sources. Neglecting the surface 
integral leads to the expression 

Ub = ^ [ B'd'r (32) 
8tt JV 

for the magnetostatic energy purely in terms of the magnetic field. 
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To test the several forms for the magnetostatic energy, we consider a 
current distribution j in a uniform magnetic field Bq. The vector potential 
for a uniform magnetic field is given by 

Ao = iBoXr. (33) 

Then the magnetostatic energy is 

f/jAo = ^/^j-(BoXr)rfV = ^Bo- l^rxjrfV = /x-Bq, (34) 

where /j, is the magnetic moment of the current distribution. 

We see that in a uniform magnetic field the magnetostatic energy depends 
only on the magnetic dipole moment. For a "point" magnetic dipole, having 
a current distribution of negligible extent, the energy in a spatially varying 
field would be given by 

f/jA = /x.B(r), (35) 

where B(r) is the magnetic field evaluated at the position of the dipole. We 
note that the sign of the magnetic dipole energy is opposite that for the 
electric dipole. That is because the current producing the magnetic dipole 
is kept constant by a constant current source. The magnetic dipole would 
rotate so as to increase its energy, thus tending to align with the magnetic 
field. For this case, with the current kept constant by a constant current 
source, the force on the point dipole will be the positive gradient of the 
energy: 

F = V(/x-B). (36) 

The magnetic dipole force thus has the same form and sign as the electric 
dipole force. 

What happens if we use Eq. (1251) to evaluate the magnetostatic energy of 
a point magnetic dipole in a uniform magnetic field? The magnetic field due 
to a point magnetic dipole is given by 

B = ^^^^'^^^3^ ~ (M-f)f] S{r). (37) 

For a magnetic dipole in a uniform magnetic field, the volume integral in Eq. 
(128D gives 

^Bo = / [a* - (/^-r)?] •Bo5(r)rfV = ^^-Bq, (38) 
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which is not the correct energy. The surface integral for this case gives 



Usurf = ^^dS-(BjXAo) 



Arc Js 

— /dS-^f^^^'^^^"^^ x(BoXr} 



Stt J 

= __L /■ dflr. [Bo(Ai-f) - f(M-Bo)] 

= ^M-Bo. (39) 

As was the case for the electrostatic energy, adding the surface integral to the 
volume integral gives the correct value /x-Bq for the magnetostatic energy. 

As was done for the electrostatic energy, the magnetostatic energy could 
also be inferred from Eq. f|T2l) to be given purely in terms of B^ by 

U^ = ^ f B^rfV, (40) 



57r Jv 

in agreement with Eq. fl32l) . This would also give the wrong answer |/x-Bo 
if used to calculate the energy of a magnetic dipole in a uniform magnetic 
field. 

We now show (as in the electrostatic case) that the deficiency in mag- 
netostatic energy if the integral of B^ is used occurs because energy would 
escape through the surface of the volume during the build up of the magnetic 
field.. We consider a uniform magnetic field Bq building up from zero to its 
final value in the presence of the magnetic dipole ^t, and calculate the energy 
that escapes through the bounding surface (assumed as a sphere of radius 
R) by integrating the Poynting vector over the surface. This gives 



t^escape = ^ ^ ^ dS-(Eo X B^) 



^^jdt<fdS-[{dtAo)xBfj,] 



-^^dS-(AoXB;x) (41) 



The vector potential Aq for the uniform field Bq is given by Eq. ( l33ll . 
and the magnetic field of the magnetic dipole by Eq. ( 1371) . Using these in 
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Eq. dUD leads to 



escape 



-— ^dS-[(BoXr)xB;x] 

f i?^(inf .[f (Bo-B^) - Bo(B;^-f )] 
_i_ ^dfii?3[(Bo-B^) - (f.Bo)(B;x.f)] 
-i-^dl]{(Bo-[3f(M-f) - m] - (Bo-f)[3(Ai-f) - (M-f)]} 
^^dfi[(Bo-/x)-(Bo-f)(/^-f)] 
^M-Bo. 



(42) 



This means that if the integral of B^ is used to determine the magnetostatic 
energy for this case, one third of the energy would escape through any bound- 
ing surface. This energy must be included, using the surface integral in Eq. 
(!27|) . to get the correct total energy. 

What combination of currents and fields can we interpret as a magneto- 
static energy density? As indicated by the steps in Eq. fl26p . the integral of 
j-A over a volume gives the magnetostatic energy within that volume. Thus 



can be considered the magnetostatic energy density. On the other hand, 
integrating B^ over a finite volume only gives the magnetostatic energy within 
the volume if the surface integral vanishes. Only in that case does the integral 
over B^ gives the same answer as that over j-A. In any event, the first 
equality in Eq. (30) shows that there is no magnetostatic energy in any region 
without current. The fairly common belief that there can be magnetostatic 
field energy in regions without current is wrong. Our conclusion is that 
B^/Stt cannot be considered a magnetostatic energy density. 

4 Magnetostatic Energy of a Permanent Mag- 
netic Dipole 

The discussion above was for a current distribution held constant in time by 
a constant current source that provided energy. For strictly classical physics. 



UjA = -j-A 



(43) 
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if magnetism is assumed to come only from such currents, this would be 
a complete description. However elementary particles (leptons and quarks) 
have permanent point magnetic dipole moments that arise from relativistic 
quantum mechanics. These permanent magnetic moments do not arise from 
electric currents, and remain constant without the influence of an external 
source. Atomic and nuclear states with quantized spin and orbital angu- 
lar momentum will also have permanent magnetic dipole moments. Then, 
a macroscopic object, such as a permanent magnet composed of polarized 
atoms locked in place, can have a permanent magnetic dipole moment. 

The force on a permanent magnetic dipole would be expected to be the 
same as for a dipole connected to a constant current source, and so be given 



This also agrees with experiment for the force on a permanent magnetic 
dipole. In fact, Eq. fHil) should be considered a separate force law for per- 
manent magnetic dipoles on a par with the Lorentz force on charge or the 
j X B force on currents, since it cannot be derived from those force laws. 

As was the case for the electric dipole force, the magnetic force on a per- 
manent magnetic dipole should equal the negative gradient of the magnetic 
energy of the dipole. Thus the energy of a permanent magnetic dipole in an 
external magnetic field Bq will be given by 



This has the opposite sign from the energy expression for a classical current 
loop but, as we have seen, that leads to the same force for each. Since the 
permanent magnetic dipole is not related to a classical current, we cannot 
derive any expression for the energy purely in terms of the fields, as in Eq. 
fl27|) . The energy of a magnetic dipole in a magnetic field is given by Eq. 
fH5|) . There is no consistent definition of energy within the magnetic dipole, 
even if it is finite in size. This is because of the inelasticity of the process of 
its magnetization. 

5 Electromagnetic Energy 

We now derive general forms for the electromagnetic (EM) energy for charge 
and current densities, and electric and magnetic fields that may be time vary- 
ing. In this section, we calculate the EM energy in the presence of polarizable 



by 




(44) 



Ufx = -Ai-Bo. 



(45) 
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material. We start with the power input to matter by an electromagnet field 
in a volume V, given by Eq. fill I) 



dUr 



matter 




f yf^^+ld^A) d\ 



dt 



Jv \ c 




The last step above assumes linear media so that 



0(atp)=p(9i0) = ^9i(p0) and y{dtA) = A.{dd) = ^dt{j-A), (47) 



which hold if and A satsfy Poisson's equation. These relations also hold for 
non-isotropic media as long as the tensors of permittivity and permeability 
are symmetric. 

The surface integral in Eq. ( 146|) will vanish for contained current distri- 
butions, provided that the surface integral encloses the entire current distri- 
bution. It will also vanish for infinitely long wires in an infinite volume if the 
scalar potential vanishes at infinite radius fast enough. Then, conservation 
of the sum of the energy in matter and electromagnetic energy implies that 
the EM energy in the volume V is given by 



In a finite volume which does not encompass all of the current, the energy 
carried out of the volume by currents in the surface integral in Eq. P6|) would 
have to be considered. 

We see that, in terms of time varying charge- current distributions and 
potentials, the expression for the EM energy is just the sum of the elec- 
trostatic and magnetostatic energies. Note that the presence of polarizable 
matter never entered the above derivation. Polarized matter does affect the 
EM energy through its affect on the scalar and vector potentials, but it does 
not affect the form of the energy density. 




(48) 
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The EM energy in a volume V can be put in terms of the EM fields using 
Maxwell's equations: 



EM 



1 

1 



svr Jv 
1 



57r Jv 



)0+ ij-A 
0(V-D) ^ 



A- (vxh 

D-V0 + V-(HxA 



c 



-H-(VxA) 



-A.(9tD) 



- (f dS-(0D + HxA) + — / (E-D + B-H)ciV 

Stt Svr Jv 

--^ / (D.^tA- A.aiD)(iV. 

Sttc Jv 



(49) 



The time derivative integral above is new. It vanishes for time constant fields, 
but the contribution of that integral and of the surface integral would have 
to be included in those cases for which they didn't vanish. In any event, a 
correct evaluation of these last three integrals would have to reproduce the 
starting point, that is Uem as given by Eq. fllSl) . The commonly used relation 



(50) 



Ueb = ^ [ (E-D + B-U)£r 

OTT Jv 



only equals the EM energy when the surface and time derivative integrals in 
Eq. (jin]) vanish. 

We look again at Eq. ([T^ . rewritten for polarizable matter as 



at 8tt Jv 



[ 9t(E-D + B-H)c/\ = -— i" dS-(ExH 

Jv An Js 



with the Poynting vector for EM energy flow given as 



An 



(ExH) 



(51) 



(52) 



Although Eq. ([5lD has been used to identify (E-D + B-H)/87r as an EM 
energy density, we see that the Poynting vector is not related directly to the 
EM energy, but only to its time derivative. If Eq. flSTj) is integrated over time 
starting with no EM energy, then some of the energy put in by matter can 
escape through the surface due to the surface integral of the Poynting vector. 
This means that, just as in the static cases, if we try to use (E-D + B-H)/87r 
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as the EM energy density, its integral may not give the total energy put into 
the volume. This is why the additional integrals in Eq. fl49l) have to be 
included to get the total energy put into the volume. 

For integrals over finite volumes our conclusion is the same as in the static 
cases. Physical EM energy can only exist within a volume if the volume 
contains charge or current. There can be no EM field energy in a volume 
that does not contain charge or current. This means that the quantity (E^ + 
B^)/87r cannot be considered an EM energy density. If EM energy resides 
anywhere, it must be within charge or current distributions. 

The point of view is sometimes taken that (E^ + B^)/87r and |(p0+j-A) 
are alternate expressions for the EM energy density, and either one can be 
chosen. This could work if both were considered as mathematical constructs 
with no direct physical significance. But, if the EM energy within a finite 
volume is an actual physical quantity, only one of the two choices can be 
applicable. A physical energy cannot have different locations, depending on 
our whim. Our preceding discussions for the static cases and the derivation 
of Eq. (H6ll show that only |(p</' + j-A) can be considered a physical EM 
energy density. 

What about the common use of the Poynting vector 7-* of Eq. (!52|) to rep- 
resent the intensity of an EM wave? To understand that, we have to enter a 
somewhat philosophical discussion about the difference between mathemat- 
ically useful quantities and physical quantities. The use over many years of 
the Poynting vector to calculate angular distributions of the emission and 
scattering of EM radiation has given correct predictions for the expected 
detection of the radiation. However, the interpretation that the radiation 
passing through the bounding surface of a volume is given differentially by 
'P-dS seems inconsistent with our discussion above. That is, we have shown 
that EM energy cannot leave a volume unless there is a charge or current 
distribution outside the volume. This means the energy cannot actually ap- 
pear outside the volume until an antenna (or retina) detects the radiation 
by means of induced currents. The Poynting vector correctly gives the prob- 
ability of radiated energy being observed in a particular direction, but the 
energy does not leave the volume until it registers in a detector. 

The situation is similar to that in quantum mechanics, where the square 
of the wave function gives the probability of a particle interacting at a par- 

•^The similarity between Classical Electromagnetism and Quantum Wave Mechanics is 
discussed in Sec. 16.4 of Ref. [1]. 
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ticular position, but the particle is not physically there until it is observed. 
In either case, the Poynting vector or the wave function (or probability cur- 
rent) is a useful mathematical construction, but it does not give the physical 
position of anything until the wave interacts with something. 

6 Conclusion 

Our conclusion, from the specific and general cases studied above, is that 
electromagnetic energy can exist only within matter, and there is no EM 
energy density in empty space. The integral of (E-D + B-H)/87r only gives 
the EM energy within a volume if the surface and time derivative integrals 
in Eq. fH^ vanish, in which case it gives the same result as the integral of 
|(p0 + j-A). This means the only valid choice for an EM energy density 
is ^(p0 + j-A). As discussed above, the Poynting vector correctly gives the 
probability of radiated energy being observed in a particular direction, but 
the energy does not leave the volume until it registers in a detector. 
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